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I. INTRODUCTION 



The operator method has been widely used as an elegant analytical tool in quantum 
mechanics jl], 0] and quantum statistics 0] for studying exactly solvable models. Since 
the publication of a remarkable paper by Delbecq and Quesne H, the operator method 
has been extended to many other possible physical applications and realizations pHnif. 
Indeed, the generalization of the notion of "shift operators" or "ladder operators" as a 
spectrum generating algebra can be effectively studied both mathematically and physically 
in terms of "nonlinear algebra" introduced by Delbecq and Quesne [Q. In fact, Chen et 
al [12| has recently reproduced the results of Bethe Ansatz |i4Hi6j for the XXX model 
through the shift operator method. A general observation in these works appears to be a 
possibility of treating the well-known technique of the algebraic Bethe Ansatz for quantum 
many-body problem as a special case of the shift operator method. Thus, all previous 
works evaluated using the algebraic Bethe Ansatz should hopefully be possible using a 
nonlinear algebraic method via shift operators. A further implication of this notion is 
the emergence of an underlying Yangian algebra. Furthermore, the factorization method 
employed in supersymmetric quantum mechanics (SQM) |HJ can also be deemed as a special 
case of the nonlinear Delbecq-Quesne algebra. 

In many realizations of physical models, the operator method can provide us with an 
alternative and a clearer picture regarding the analytical determination of energy eigenstates 
and eigenvalues. Specifically, it has been found to be a useful tool in the solution of bound 
state problems such as the string postulate in Bethe Ansatz or the raising and lowering of 
energy levels in SQM. As in Bethe Ansatz and SQM, knowledge of a certain reference state, 
like the highest weight state in Bethe Ansatz or the ground state ansatz for the superpartner 
potential in SQM, in the operator method |l] permits the analysis and determination of the 
full spectrum of a physical system. Moreover, as mentioned in ref . |2| , the knowledge of the 
ground state wave function also determines the potential and hence the Hamiltonian of a 
system up to a constant value. However, in principle given the Hamiltonian of a physical 
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system, surely it should be possible to determine all the energy levels. In this paper, we 
show that provided certain consistency relations are satisfied, the entire spectrum of a given 
Hamiltonian can be determined by shift operator method even if there is no prior knowledge 
of a reference state. 

To illustrate the power of the shift operator method, we apply the technique to one- 
dimensional oscillator with anharmonic potentials. The study of anharmonic potential 
||18|^0| has always been an exciting and interesting field due to its broad applications in 



quantum field theory |21|,|22||, nuclear models, atomic and molecular physics |23 H25fl , con 



densed matter physics |26],[27]], statistical physics and chemical physics [p8fl . Indeed, nu- 



merous numerical methods including renormalized strong coupling expansion, perturbation 
expansion, supersymmetric quantum mechanics, WKB, iteration based on the generalized 
Bloch equation, state-dependent diagonalization, Hill determinant method, phase-integral 
approach, iterative Bogoliubov transformations, eigenvalue moment method, perturbative- 
variation, and algebraic method have been proposed to investigate these anharmonic po- 
tentials pU|-|^5[]. However, due to its inherent intractability, no analytic solution for energy 



spectrum and eigenstates has been obtained so far except in some special cases. Indeed, for 
the first time, we have partially resolved some of these technical issues and obtained analytic 
solutions for a wide class of anharmonic potentials under certain conditions. 

It is instructive to see how SQM can be regarded as a special case of the nonlinear 
algebraic operator method. Before we proceed further to show this point, we first recall 
some relevant definitions in operator method and nonlinear algebra |jl|,|4, LT|,|13|| . Let H be 



an observable satisfying the eigenfunction equation H\ip >= E\tp >. The operators L ± 
satisfying 

[H,L ± ]=L ± f ± (H), (1) 

are called "shift operators" of H. In Eq. ([!]), f ± (H) are real functions of H so that the 
values of its action on eigenstates of H are reals. Thus, these values can be interpreted as 
energy gaps if H is the Hamiltonian of the system. Nevertheless, it is not necessary for the 



mutually adjoint condition 

(L+)t = L~ 

to hold and there is generally no constraint on the commutation relation between L + and 
L~, i.e. [L + ,L~]. However, if the mutually adjoint condition is satisfied, we obtained the 
"nonlinear algebra" |4| defined by the relations 

[H,L+] = f(H)L + , 
[H,L-] = -L~f{H), 
[L-,L + ]=g(H). 

A simple rearrangement of the above relations gives 

{H-f{H))L + = L + H, 

L-(H-f(H)) = HL-. (2) 



It is interesting to compare Eq.(||) with the analogous relations in SQM |T7 



HtQ + = Q + H 2 
Q-H x = H 2 Q- 

Q- = (Q + ) t (3) 



and identify 



Q + ~ L + Q- ~ L~ 
H X ~H - f(H) H 2 ~ H. (4) 

Hence, one can always regard SQM as a specific type of nonlinear algebraic method. 

The purpose of this paper is to describe the operator method fl1,|TT|-|T3|| and apply it 



to study the one-dimensional oscillator with anharmonic potentials. In section [TT], we show 
how a shift operator method can be formulated using a simple instructive example of the 
harmonic oscillator. The general procedures to solve anharmonic potentials are presented in 



section [II. In section [IV], we solve the relevant equations for the one- dimensional oscillator 
with cubic-quartic potential and show that we can get the exact analytical expression for the 
ground state energy provided certain consistency relations are satisfied. We then compare 
our analytical result with some numerical results in section |V|. Our analytical expressions 
can be used to verify the accuracy of existing numerical results. We conclude with some 



brief remarks in section VI 



II. APPLYING OPERATOR METHOD TO THE HARMONIC OSCILLATOR 



In this section, we apply the operator method to the harmonic oscillator. The Hamilto- 
nian of the harmonic oscillator is given by 

1 



H=-(p 2 + x 2 ). 



(5) 



Together with the well-known commutation relation between x and p, i.e. [x,p] = i (h = 1) 
we get 



Using the notation 



[H, x] = —ip 
[H,p] = ix. 

[H, (x,p)} :=([H,x],[H,p]), 



(6) 



Eq.(O) can also be rewritten as 



( A 



[H, (x,p)) = (x,p) 



i 



{x,p)M 



(7) 
(8) 



where, the matrix, M, is called the "coefficient matrix". To diagonalize the coefficient 
matrix, we look for a transformation, U, such that M = UDU^ 1 where D is a diagonal 
matrix. In general, provided all the entries of U commute with Hamiltonian H, the entries 
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of the matrix U are functions of c numbers, the Hamiltonian or any conserved quantity I of 
the system. Thus, we have 



[H, (x,p)U(c, H, /)] = (x,p)U(c, H, i)U-\c, H, I) 

( 



d A 



\~ i0 J 



(x,p)U{c,H,I) 



\ x {c,HJ) 
\ 2 {c,H,I) 



U(c,H,I) 



(9) 



It is easy to see that 



U = -r 



i A 



so that Eq.(^) can be rewritten as 



[H, (X'P)^ 



v 



-i % 



i A 



v -. z y 



1 
-1 



(10) 



Identifying 



[a , a) 



id A 



\- 1 1 ) 



we obtain 



[H,a+] = a + 
[H, a] = —a. 

The "creation" and the "annihilation" operators are thus obtained naturally. Following the 
literature, we can also call them "shift" or "ladder" operators. Although, the application 
of the shift operator to harmonic oscillator is very simple, it can still serve an illustrative 
example. 

For the harmonic oscillator, a closed algebra is obtained. This closure in the algebra 
enables the whole spectrum of the system to be determined using only one pair of shift 
operators, in which one of the operator raises while the other lowers the energy levels. In 
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general, the existence of a closed algebra is not assured, and so, we usually have to deal with 
an unclosed algebra with infinitely many shift operators. Thus, we have to use infinitely 
many shift operators pairs to generate the whole spectrum. Moreover, corresponding to each 
energy level, we have a pair of distinct shift operators in which one raises while the other 
lowers the energy level. 

It is interesting to note that, the eigenvalues of the coefficient matrix correspond to the 
amount of shifted energy associated with the various shift operators. And, except for some 
special Hamiltonians, the coefficient matrix, its eigenvalues and the transformation matrix, 
are all functions of the Hamiltonian. In nonlinear algebraic method, the energy shifts are 
no longer uniform and the gap between any two adjacent energy levels is related to their 
relative positions in the spectrum. Thus, we find that shift operator method can provide us 
with a clearer physical picture for the mapping of the energy level in a physical system. 

Note that we can only extract the energy gaps rather than the energy levels. Moreover, 
we have no information regarding to the energy of the ground state |Ij]. However, in some 
cases, it is possible to determine the exact energy of the ground state. This last result 
constitutes the gist of one of the most important aspect of our paper. 

III. OPERATOR METHOD TO ANHARMONIC OSCILLATOR: THE GENERAL 

PROCEDURE 

In this section, we confine ourselves to the one- dimensional harmonic oscillator in which 
the anharmonic potential contains cubic and quartic terms. Nevertheless, the procedure 
developed here is very general and can be applied to other types of anharmonic potentials. 
The Hamiltonian is given by 



H 



dx 2 



+ ax 2 + /3x 3 + 7X 4 , 



(11) 



and it is easy to show that 




(12) 
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[H, x n i = -n(n - l)x n - 2 4- + 2nx n - l H 
dx dx 

-2a(n + l)x n+1 - 2p{n + ^)x n+2 



2j(n + 2)x 



n+3 



Rewriting Eqs. (|l2| , |T3D into matrix form, we have 



^-/^/"^ {^X ^ OC i ^ ^ " " " « . */] — } *^ 1 *^ 5***5 7 " / 



2 a3 



. d d 2 d d \ -\t 

U/iL> LLJU KiitXj U/iL> 



L" > V. J > ^ 7 ' ^ 7 ' 

ax ax ax 



,x" 



d 

dx' 



f /y ry ry ... T* 



■)M 2 



, d d 2 d n d NAT 

ax dx dx dx 



where 



M 1 



0-60 
-12 
-20 





-n(n — 1) 



-2 
0-4 
0-6 





-2n 



L 1 = (0, -2, 0, 0, 0, •••) 



S 



and 



M, 



-2a 





AH 










-3/3 


-4a 





AH 







_4 7 


-5(3 


—6a 





8H 







—67 


-1(3 


-8a 


"■ 










-87 


-9/5 


-10a '■ 


2nH 


n 


n 


n 

u 


— 1 Do/ 


— 11/9 ' • 


n 

u 














-127 "• 


. -2(n + l)a 














• 


• -2(n + l){3 














"■ 


■ -2(n + 2) 7 



(19) 



No 



0-20 
0-6 
-12 





-n(n + 1) 



(20) 



L 2 = (0, 2H, 0, 0, 0, •••) 



(21) 



It is interesting to note that for other types of potentials, all matrices except for M 2 are 
exactly the same. 

From Eqs.(|14,15), we find 



\fj (•r <r 2 <f 3 . . . <f n ...Iffi. (A. r l A. ... d_ \ q] 

[£1,{X,X ,X , ,X , / rl ~r V ' cix ' ' ' X das' 

= (x, x 2 , x 3 ,---,x n r - ■)(M 1 R + M 2 S) + (£, z^, x 2 -£, ■ ■ • , ■ • ■) (7V 1jR + AT 2 S) + (^i? + L 2 S) 



(22) 



Thus, if the following relations exist 



M ± R + M 2 S = RT 
N X R + N 2 S = ST 



(23) 
(24) 



then we have 



dx ' 



■)s\ 



(25) 



\fj f A A2 A3 , , , An ~\ P I ( d d ™2_rf_ ... „r, 

= {(x, x 2 , £ 3 , ■ • • , x n , ■ ■ -)R + {£, x£ x 2 ±, ■■■,x n £-- -)S}T + {L X R + L 2 S) 

To diagonalize the coefficient matrix T, let us suppose that C/ is the transformation matrix 
needed. The eigenvalues and the corresponding shift operators can therefore be written as 



A = U~ L TU = diag(Ai, A 2 , A 3) • • • , A n , • • ■) 



(26) 



and 



(A 1 ,A 2 ,A 3 ,---,A n ,---) 

= {{x, x\ x\---,x n ,-- -)R + (£ x£ x 2 ±, • ■ • , x n ± ■ ■ -)S}U + (LiR + L 2 S)UA~ 1 

respectively. In fact, Eq.(|25D can be written into a more succinct form as 

/ 



(27) 



[H, {A 1 ,A 2 , A 3 , ■ ■ ■ , A n , ■ ■ •)] = {At, A 2 , A 3 , ■ ■ ■ , A n , 



Ai 
A 2 
A 3 







Ar. 



v • 



, (28) 



: / 



or 



[H,A 1 ] = AiAi 
[H,A 2 ] = A 2 X 2 



[H, A n ] = A n X 
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Thus, in general we get an unclosed algebra with infinitely many shift operators. In case 
of the usual harmonic oscillator, we have a pair of shift operators which can generate the 
entire spectrum. In the present case, we have infinitely many pairs of shift operators in 
which each pair is responsible for raising and lowering the corresponding energy level. The 
whole spectrum can only be generated by the infinite set of the shift operators acting on the 
ground state. Hence, it is natural to see that all the eigenvalues of the coefficient matrix 
are only dependent on the energy of ground state. That is, Aj, (i = 1, 2, 3, • • • , n, ■ ■ •) are 
the functions of the ground state energy. Therefore, we can identify the operator H in the 
coefficient matrix T as the ground state energy! Based on this observation, we can then 
get the analytic expression for the energy of the ground state provided certain consistency 
relations hold. 

It remains to solve for the matrices R, S and T given the matrices M\, M 2 , Ni and N%. 
Once we have obtained the matrices R, S, T and diagonalize T, we effectively obtained all 
the shift operators and we can then reconstruct the full spectrum. It is instructive to note 
that one of the matrices R or S is redundant and we set it to unity. For convenience, let us 
set S in Eqs. (|2"3"|,|2~3D to unity so that the problem is reduce to the solution of R, T from the 



following equations 

M 1 R + M 2 = RT (29) 
N X R + N 2 = T (30) 

and the diagonalization of T. Note that although we have considered a specific form for 
the potential, the procedure developed here is very general and can be applied to more 
complicated cases. 

IV. ANALYTIC EXPRESSION FOR THE GROUND STATE ENERGY: A CLASS 

OF ANHARMONIC POTENTIALS 

To solve Eq.(|2lf) and Eq. (|3"0|) , we first make the following observation. If we define 
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G 



1 
2 
3 







n 



100 

0010 

0001 

■•• 1 ••• 



Q 



■J 



1 
1 
1 





(31) 



(32) 



(33) 



we have 

Mi = -PGPG 

M 2 = -2aG + 2HPGP -2-iQGQ- !3{2G- l)Q 
Ni = -2G 
N 2 = -GPGP. 

Here, the matrices G, P and Q are nothing but representations of particle number, creation 
and annihilation operators for the usual harmonic oscillator in Fock space. It is easy to 
check that 



[G,P] = -P, [G,Q}=Q, PQ = l, QP = diag(0,l, !,•••)• 



(34) 



Furthermore, if we set 



then Eq.pO) becomes 



W = NiR, 



T = W + N 2 R = N^W. 
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From Eq.(||), we notice that N^N^ 1 = N 2 , and as a consequence, we can combine 
Eq.(p9D and the above expression for the matrix T, to yield the expression 

W 2 + [W, N 2 ] - A^!M 2 = 0. (35) 

Unfortunately, due to the infinite dimensionality of the matrices, the solution, except in 
some special cases, is not known. 

In the case of the one-dimensional oscillator with cubic-quartic term, if we let W assume 
the following form 

W = 2(aG + bGQ + cGP), 
it is not difficult to show using Eq.(p^) that 

W 2 + [W, N 2 ] - NM = 4{(a 2 + 2bc - a)G 2 + (06 - \(3){2G 2 - G)Q + (b 2 - >y)GQGQ 

{ac + b){2G 2 + G)P + (c 2 + a + H)GPGP} (36) 

so that consistency naturally requires the following conditions 

a 2 + 2bc - a = (37) 

ab - -p = (38) 

b 2 - 7 = (39) 

ac + b = (40) 

c 2 + a + H = (41) 



in order that Eq.fl3"5D be satisfied. In particular, Eq.(f4"TD gives the ground state energy of 
the system. 

As explained in the previous section, we have equated H to one of its eigenvalues, namely 
the ground state energy. Thus, in order that the full spectrum be generated from the in- 
finitely many raising operators acting on the ground state, so that each raising operator 
generates its own corresponding energy level via its action on the ground state, the consis- 
tency relations in Eq.(^)-Eq.(^) must hold. 
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It is easy to see that the solution to Eq.(|3~7l)-Eq.(f41~|) is 



' ! Vl c=-2l (42) 



2^ v ' P 



H = - JL (43) 



and 



P 2 2^7 



B 2 -y 3 / 2 

7- " 4 V " « = ( 44 ) 



' b=-^ c =-2^ ( 45 ) 



2^ v /3 

B 2 7 3 / 2 

|- + 42 r -a = 0. (47) 

The additional constraint Eq.(|44|) (or Eq.([|7])) seems to indicate that the ansatz for W 
may be too simplistic. However, we do not have a direct solution of W from Eq.(|35[) for 
the given N\, N2 and M2. Despite all these, we can still obtain an analytic result for exact 
energy of the ground state for the given anharmonic potential and this result can prove to be 
valuable for analyzing the accuracy of existing numerical methods. Besides, as mentioned 
before, it is anticipated that the operator method can provide us with an analytic tool for 
investigating the ground state PJ], something which is not possible using existing numerical 
approaches. 

To compute the energy levels of the excited states, we need to diagonalize the following 
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infinite dimensional matrix 



/ 

2a 


2c 


-2 












46 


4a 


4c 


-6 












66 


6a 


6 c 


-12 












86 


8a 


8c "■ 















106 


10a "■ 


• —n{n + 1) ■ 
















126 '• 


2nc 
















'■ 


2na 



















2n6 




V ; 















with a, 6 and c are given by Eq . (f42|) or Eq.([45[). Once we have diagonalized this matrix, 
using the result of the ground state energy, all the other energy levels can be obtained. 
Unfortunately, we do not know any method for diagonalizing it at present moment despite 
the apparent simplicity and symmetry in the matrix. 

V. COMPARISON OF NUMERICAL WITH THE ANALYTIC RESULTS FOR 

GROUND STATE ENERGY 

In this section, we use our analytic results to check against previous numerical com- 
putations. For numerical simulation, we use the method proposed by Ho et al. using the 
state-dependent diagonalization method. As claimed in ref. ]30]| , this method is very accu- 
rate and efficient compare to other numerical methods for calculating the energy eigenvalues 
and eigenfunctions of the one-dimensional harmonic oscillator with anharmonic potentials. 

Due to the consistency requirements in Eq . (f44"D or Eq.(f47|), the three coefficients a, j3 and 
7 in 

V(x) = ax 2 + (3x^ + 7a; 4 

are not mutually independent. It is not difficult to solve for (3 using Eq.(JP]) or Eq.(f4"T|) for 
given values of a and 7. The final solution is 
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a 
3 



iH 7 + V 72 - ( f )3 



(48) 



corresponding to Eq. (fPj) and —/3 to Eq.(f4"?|) with the ground state energy being same in 
both cases. Specifically, the ground state energy is given by 



The dependence of (3 on a and 7 is shown in Fig.(|l]) 



(49) 



Setting a — 1, —2, the values of E as a function of 7 for the analytic formula and the 
numerical simulation using the state dependent diagonalization |30| are shown respectively in 
Fig.(|) and Fig.(|3|)) with the corresponding potential for particular values of 7 as a function 
of x shown in Fig.(|4])). The difference in the energy computed from the numerical simulation 
and the analytic formula as a function of 7 is also plotted in Fig.([|) and Fig.(^) for a = 1 
and a = —2 respectively. As shown in the figures, the numerical simulation using state- 
dependent diagonalization is in excellent agreement with our analytical formula. Moreover, 
we see that the state-dependent-diagonalization method provides sufficiently high accuracy 
for all practical purposes. 



VI. CONCLUSIONS 

In summary, we have formulated a general procedure for extracting the full spectrum of 
a physical system with arbitrary potentials using the operator method. We have applied it 
to solve the one-dimensional harmonic oscillator with anharmonic problem. The analytic 
expression for the ground state energy is obtained for a large class of anharmonic potentials. 
Our results can be used to verify the accuracy of existing numerical methods. However, 
in order to get the full spectrum, we need to diagonalize an infinite dimensional matrix. 
This last problem is to our knowledge an open mathematical question except for some very 
special cases. 
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Our method has also confirmed the notion that once the Hamiltonian of certain system 
is given, all the energy levels can be determined through the operator method even without 
any prior knowledge of a reference state, provided certain consistency relations hold, a feat 
considered impossible hitherto. 

Finally, from Eq.QSSD and Eq.fl4~9|), it easy to see that in the limit when 7 — > 00, E — > 
— 1(27) 1//3 , confirming earlier analysis on the mathematical properties of ground-state energy 



obtained in ref. [36. 46 
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FIG. 1. The dependence of [3 on 7 and a. 



FIG. 2. The comparison of our analytic result (bold line) with that obtained by state-dependent 
diagonalization method (represented by the data points) for a = 1. 

FIG. 3. The comparison of our analytic result (bold line) with that obtained by the 
state-dependent diagonalization method (represented by data points) for a = —2. 

FIG. 4. The potential, V(x), as a function of x for (a) a = 1, 7 = 1/2 and (b) a = —2, 7 = 1/4. 

FIG. 5. The difference in the energy between the analytic and the numerical data as a function 
of 7 for a = 1. 

FIG. 6. The difference in the energy between the analytic and the numerical data as a function 
of 7 for a = —2. 
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